A new hybrid LDA and Generalized Tight-Binding method for the electronic 
structure calculations of strongly correlated electron systems 
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A novel hybrid scheme is proposed. The ab initio LDA calculation is used to construct the Wannier 
functions and obtain single electron and Coulomb parameters of the multiband Hubbard-type model. 
In strong correlation regime the electronic structure within multiband Hubbard model is calculated 
by the Generalized Tight-Binding (GTB) method, that combines the exact diagonalization of the 
model Hamiltonian for a small cluster (unit cell) with perturbation treatment of the intercluster 
hopping and interactions. For undoped La2Cu04 and Nd2Cu04 this scheme results in charge 
transfer insulators with correct values of gaps and dispersions of bands in agreement to the ARPES 
data. 
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I. INTRODUCTION 

A conventional band theory is based on the density 
functional theory (DFT)i and on the Local Density Ap- 
proximation (LDA)— within DFT. In spite of great suc- 
cess of the LDA for conventional metallic systems it ap- 
pears to be inadequate for strongly correlated electron 
systems (SCES). For instance, LDA predicts La2Cu04 
to be a metal whereas, in reality, it is an insulator. Sev- 
eral approaches to include strong correlations in the LDA 
method are known, for example LDA+Lfci and LDA- 
SIG— . Both methods result in the correct antiferromag- 
nctic insulator ground state for La2Cu04 contrary to 
LDA, but the origin of the insulating gap is not correct. 
It is formed by the local single-electron states splittcd 
by spin or orbital polarization. In these approaches the 
paramagnetic phase of the undoped La2Cu04 (above the 
Neel temperature Tm ) will be metallic in spite of strong 
correlation regime U S> W, where U is the Hubbard 
Coulomb parameter— and VP is a free electron bandwidth. 
The spectral weight redistribution between Hubbard sub- 
bands is very important effect in SCES that is related to 
the formation of the Mott-Hubbard gap in the param- 
agnetic phase. This effect is incorporated in the hybrid 
LDA+dynamical mean field theory (DMFT) (for review 
see RefAL^) and LDA++ approaches — The electron self- 
energy in LDA+DMFT approach is calculated by the 
DMFT theory in the limit of infinite dimensiorji2iii and 
is k- independent, E k (E) -> E (£ , )i2J^ That is why the 
correct band dispersion and the ARPES data for High- 
T c compounds cannot be obtained within LDA+DMFT 
theory. Recent development of the LDA+cluster DMFT 
methodi^ii^ and spectral density functional theory^ 
gives some hopes that non-local corrections may be in- 
cluded in this scheme. 

A generalized tight-binding (GTB)i£ method has been 
proposed to study the electronic structure of SCES as a 
generalization of Hubbard ideas for the realistic multi- 



band Hubbard- like models. The GTB method combines 
the exact diagonalization of the intracell part of the 
Hamiltonian, construction of the Hubbard operators on 
the basis of the exact intracell multielectron eigenstates, 
and the perturbation treatment of the intercell hoppings 
and interactions. A similar approach to the 3-band p — d 
model of cupratesiSiiS is known as the cell perturba- 
tion methodi 20 i 21 ' 22 The practical realization of the GTB 
method for cuprates required an explicit construction of 
the Wannier functions to overcome the nonorthogonal- 
ity of the oxygen molecular orbitals at the neighboring 
CuOq cells— ^ The GTB calculations for undoped and un- 
derdoped cuprates are in good agreement with ARPES 
data both in the dispersion of the valence band and in the 
spectral intensityi^&Si A strong redistribution of spectral 
weight with hole doping and the formation of the in-gap 
states have been obtained in these calculations. Simi- 
lar GTB calculations for the manganites has been done 
recently—^ 

As any model Hamiltonian approach the GTB method 
is not ab initio, there are many Hamiltonian parameters 
like intraatomic energy levels of p and d electrons, various 
p — d and p — p hopping parameters, Coulomb and ex- 
change interaction parameters. These parameters have 
been obtained by fitting the set of optical, magnctio2£ 
and ARPES^ data. Generally the question arises how 
unique the set of parameters is. To overcome this restric- 
tion we have proposed in this paper a novel LDA+GTB 
scheme that allows to calculate the GTB parameters by 
the ab initio LDA approach. 

The paper is organized as follows: In Section [H] the 
construction of Wannier functions from self-consistent 
LDA eigenfunctions as well as ab initio parameters of 
the multiband p — d model for La2Cu04 and Nd2Cu04 
are given. A brief description of the GTB method is 
done in Section [TTU Section |IV| contains the LDA+GTB 
band structure calculations for La2Cu04 and Nd2Cu04. 
The effective low-energy t — J* model with ab initio pa- 
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rameters is presented in Section Section IVII is the 
conclusion. 



II. CALCULATION OF AB INITIO 
PARAMETERS FROM LDA 

To obtain hopping integrals for different sets of bands 
included in consideration we apply projection procedure 
using Wannier functions (WFs) formalism^ WFs were 
first introduced in 1937 by Wannier— as Fourier trans- 
formation of Bloch states |*0»ic) 



\W 



-ikT 



IV'ik), 



(1) 



where T is lattice translation vector. N is the num- 
ber of discrete k points in the first Brillouin zone and 
i is band index. One major reason why the WFs have 
seen little practical use in solid-state applications is their 
nonuniqueness since for a certain set of bands any or- 
thogonal linear combination of Bloch functions \ipik) can 
be used in 0J. Therefore to define them one needs 
an additional constraint. Among others Marzari and 
Vanderbilt— proposed the condition of maximum local- 
ization for WFs, resulting in a variational procedure. To 
get a good initial guess authors of Refl28l proposed to 
choose a set of localized trial functions \4> n ) and project 
them onto the Bloch states |^tk)- It was found that this 
starting guess is usually quite good. This fact later led to 
the simplified calculating scheme^ where the variational 
procedure was abandoned as in present work and the re- 
sult of the aforementioned projection was considered as 
the final step. 



A. Wannier function formalism 

To construct the WFs one should to define a set of trial 
orbitals \(f> n ) and choose the Bloch functions of interest by 
band indexes (Ni, . . . , N2) or by energy interval (EijE^). 
Non-orthogonalizcd WFs in reciprocal |W„k) space are 
then the projection of the set of site-centered atomic- 
like trial orbitals \4> n ) on the Bloch functions IV'ik) of the 
chosen bands: 

l^nk) = l^k)(^k|0„), (2) 

l (B 1 < £l (k)<£; 2 ) 

where 6j (k) is the band dispersion of i-th band obtained 
from self-consistent ab initio LDA calculation. In present 
work we use LMT-orbitals2fi as trial functions. The Bloch 
functions in LMTO basis are defined as 



in 1 



(3) 



where [i is the combined index representing qlm (q is the 
atomic number in the unit cell, I and m are orbital and 



magnetic quantum numbers), <rf^(r) are the Bloch sums 



of the basis orbitals ip^ (r — T) 

#W = 7=E elkT <M r - T )' 



and the coefficients are 



c % = (^vlV'ik)- 



(4) 



(5) 



Since in present work \<f> n ) is an orthogonal LMTO 
basis set orbital (in other words n in \<f> n ) corresponds 
to the particular qlm combination), then (ipik\<fi n ) = c k *. 
Hence 

N 2 N 2 

\w nk ) = £ |^ k }4* = £ E c >«*l$>- ( 6 ) 

i=Ni i=Ni A 1 

In order to orthonormalizc the WFs © one needs to 
calculate the overlap matrix O nn ' (k) 

N 2 

O nn ,(k) = (W nk \W n , k ) = £ c* c&, (7) 

i=Ni 

then its inverse square root S nn ' (k) is defined as 



s nn ,(k) = o- n V 2 (k). 



(8) 



In the derivation of (0 the orthogonality of Bloch states 
(V^klVVk) = <W was used. 

From © and (jHJ , the orthonormalized WFs in k-space 
\W n k) can be obtained as 



N 2 



|^„k) = 5>„ n ,(k)|WV k ) = J2 l^ k > 5 ™ ' 



f=JVi 



4* = (^klWnk) = ^5 nn ,(k)c^. 

n' 

Then the matrix element of the Hamiltonian H WF in 
reciprocal space is 



N 2 



k' i=Ni 



H%rw = (w nk |(^EE ^^(ko^k'iji^k) 

(9) 



— c ni c n'i e iO i )- 

i=JVi 

Hamiltonian matrix element in real space is 

, N 2 



ikT 



k i=Ni 



here atom n' is shifted from its position in the primary 
unit cell by a translation vector T. For more detailed 
description of this procedure see— 
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FIG. 1: (Color online) Comparison of the band structure 
of La2Cu04 from LDA calculation (dotted lines) and from 
projection on the Cu-d x 2_ y 2 and O p -p x , O p -p y set of orbitals 
(bold solid lines). Fermi level coresponds to zero energy. 



FIG. 3: (Color online) Comparison of the band structure 
of Nd2Cu04 from LDA calculation (dotted lines) and from 
projection on the Cu-d x 2_ y 2 , 01-p x and 01-p y set of orbitals 
(bold solid lines). Fermi level corresponds to zero energy. 
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FIG. 2: (Color online) The same as in Fig.Qbut projection is 
done on the Cu-d x 2 _ y 2 , Cu-d 3z 2_ r 2 and O p -p x , O p -p y , O a ~Pz 
set of orbitals. 



FIG. 4: (Color online) The same as in Fig.|2]but projection is 
done on the Cu-d x 2_ y 2 , Cu-d 3z 2_ r 2 and 01-p x , 01-p y , 02-p z 
set of orbitals. 



B. LDA band structure, hopping and Coulomb 
parameters for p- and n-type cuprates 

Basically all cuprates have one or more CuC>2 planes 
in their structure, which arc separated by layers of other 
elements (Ba, Nd, La, ...). They provide the carriers 
in Cu02 plane and according to the type of carriers all 
cuprates can be divided into two classes: p-type and n- 
typc. In present paper we deal with the simplest rep- 
resentatives of this two classes: La2- a: Sr a ;Cu04 (LSCO) 
and Nd2- a; Ce a ;Cu04 (NCCO) correspondingly. 

LDA band calculation for La2Cu04 and Nd2CuC>4 
was done within LMTO method^ using atomic sphere 
approximation in tight-binding approach^ (TB-LMTO- 
ASA). In the case of Nd 2 Cu04 Nd-4/ states were treated 
as pseudocore states. 



La2Cu04 at the low temperature and zero doping has 
the orthorhombic structure (LTO) with the space group 
Bmabm^ The lattice parameters and atomic cordinates at 
10 K were taken from Ref|H to be a=5.3346, b=5.4148 
and c=13.1172 A, La (0, -0.0083, 0.3616), Cu (0, 0, 0), 
O p (0.25, 0.25, -0.0084), O a (0, 0.0404, 0.1837). Here 
and below O p denotes in-plane oxygen ions and O a - 
apical oxygen ions. In comparison with high temerature 
tcragonal structure (HTT) orthorhombic La2Cu04 have 
two formula units per unit cell and the CuOg octahedra 
are rotated cooperatively about the [110] axis. As a result 
Op ions arc slightly moved off the Cu plane and four in- 
plane La-O a bond lengths are unequal. 

Nd2Cu04 at the room temperature and zero dop- 
ing has the tetragonal structure with the space group 
IA/mmmSl also called T'-structure. The lattice paramc- 
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tcrs are a=b=3. 94362, c=12.1584 KM Cu ions at the 2a 
site (0, 0, 0) are surrounded by four oxygen ions 01 which 
occupy 4c position (0, 1/2, 0). The Nd at the 4e site (0, 
0, 0.35112) have eight nearest oxygen ions neighbors 02 
at 4d position (0, 1/2, 1/4)i24 One can imagine body- 
centered T'-structure as the HTT structure of La2Cu04 
but with two oxygen atoms moved from apices of each oc- 
tahedron to the face of the cell at the midpoints between 
two oxygen atoms on the neighbouring Cu02 planes. In 
other words Nd2Cu04 in T'-structure has no apical oxy- 
gens around Cu ion. 

The LDA band structure of both compounds along 
the high-symmetry lines in the Brillouin zone is shown 
in Figs. 1 1141 by dotted lines. The coordinates of high- 
symmetry points in BZ are given on top of each picture. 
The complex of bands in the energy range (-8, 2.5) eV 
consists primarily of Cu-3d and 0-2p states. The total 
bandwidths amount 10 eV for La-cuprate and 7 eV for 
Nd-cuprate. Contribution of Cu-3d and 0-2p orbitals to 
the different bands is displayed by arrows. 

One can see that the band crossing ~Ef have character 
of Cu-d x 2_ y 2 and O p -p XiV for La2Cu04 and C\i-d x a_ y a, 
01-p x .y in the case of Nd2Cu04. It corresponds to an- 
tibonding pda orbital. So for hoppings calculation the 
projection on Cu-dj^^, O p -p x , O p -p y orbitals for La- 
cuprate and Cu-(i x 2_ y 2, Ol-p^, 01-p y orbitals for Nd- 
cuprate was done. Such set of orbitals corresponds to 
the 3-band p — d model. The bands obtained by the de- 
scribed in Sec. Ill Al projection procedure are shown by 
solid lines in Figs. Q]and|3J It is clearly seen that in case 
of La2Cu04 3-band model did not reproduce the band 
crossing Ej? properly (Fig.Q] SR direction). 

Since 3-band p — d model didn't provide proper de- 
scription of the LDA bands around Fermi level the pro- 
jection on more complex set of trail orbitals for both 
compounds was done. The resulting bands are plotted 
by solid lines in Figs. and Corresponding multi- 
band p — d model contains Cu-d x 2 _ y 2 , Cu-d 3z 2_ r 2, O p - 
Pxi Op-py, O a -p z states for La2Cu04 and Cu-d x 2_ y 2, Cu- 
d 3z 2_ r 2, Ol-p x , Ol-p y , 02-p z states for Nd 2 Cu04. The 
energy range for projection was (-8.4, 2.5) eV and (-8, 2) 
eV for the case of La-cuprate and Nd-cuprate correspond- 
ingly. The main effect of taking into account Cu-d3 Z 2_ r 2 
and O a -p z states for La2Cu04 is the proper description of 
the band structure (in comparison with LDA calculation) 
at the energies up to 2 eV below Fermi level. From Fig. [31 
and |H one can see that in case of Nd2Cu04 both sets of 
trial orbitals properly describe the LDA band crossing 
the Fermi level which has Cu-d x 2 _ y 2 symmetry. At the 
same time its bonding part does not agree well with the 
LDA bands since projection did not include all Cu-d and 
O-p orbitals. 

The resulting hopping parameters and energy of par- 
ticular orbitals for two sets of trial orbitals are presented 
in Tables [I]and[n] The second column contains the con- 
necting vector T between two sites. It is clearly seen that 
hoppings decay quite rapidly with distance between ions. 

For the multiband p — d model the values of 



TABLE I: Hopping parameters and single electron energies for 
orthorhombic La2Cu04 obtained in WF projection procedure 
for different sets of triai orbitais (aif vaiues in eV). Here x 2 , 
z 2 , Px, Py, Pz denote Cu-d x 2_ y 2, Cu-d 3z 2_ r 2, O p -p x , O p -p y , 
Oa-Pz orbitals correspondingly. The 3-d and 4-th columns 
correspond to bases of the 3-band and the multiband p — d 
models respectively. 



Hopping 


Connecting 


Cu-x 2 


Cu-a; 2 , z 2 




vector 


O Px, Py 


O-Px, Py, Pz 






E x 2 =-1.849 


E x 2 =-1.849 






E Pt =-2.767 


E z 2 =-2.074 






E Ph =-2.767 


E Px =-2.806 
E Py =-2.806 
E Pi =-1.676 


t(x 2 ,x 2 ) 


(-0.493,-0.5) 


-0.188 


-0.188 


t'(x 2 ,x 2 ) 


(-0.985, 0.0) 


0.001 


0.002 


t(z 2 ,z 2 ) 


(-0.493,-0.5) 




0.054 


t'(z 2 ,z 2 ) 


(-0.985, 0.0) 




-0.001 


t(x 2 ,Px) 


(0.246,0.25,-0.02) 


1.357 


1.355 


t'(x 2 ,Px) 


(-0.739,0.25,-0.02) 


-0.022 


-0.020 


t(z 2 ,Pi) 


(0.246,0.25,-0.02) 




-0.556 


t'(z 2 ,Px) 


(-0.739,0.25,-0.02) 




-0.028 


t(^,Pz) 


(0,0.04,0.445) 




0.773 


t'(* 2 ,Pz) 


(-0.493,-0.46,-0.445) 




-0.011 


t(Ps,Py) 


(0.493, 0.0) 


-0.841 


-0.858 


t'(Px,Py) 


(0,0.5,0.041) 


0.775 


0.793 


t"(px,Py) 


(0.985,0.5,0.041) 


-0.001 


-0.001 


t(Px,Pz) 


(-0.246,-0.21,0.465) 




-0.391 


t'(j>x ,Pz) 


(0.246,0.29,-0.425) 




-0.377 


t"(j>x ,Pz) 


(0.246,-0.21,-0.746) 




0.018 



Coulomb parameters arc also required. For Cu in 
La2Cu04 they were obtained in constrained LDA super- 
cell calculations^ to be U = 10 eV and J = 1 eV& 
For the Nd2Cu04 we will use the same values of these 
parameters. 

III. GTB METHOD OVERVIEW 

As the starting model that reflects chemical structure 
of the cuprates it is convenient to use the 3-band p — d 
modeliSiiS or the multiband p — d models While the 
first one is simplier it lacks for some significant features, 
namely importance of d z 2 orbitals on copper and p z or- 
bitals on apical oxygen. Non-zero occupancy of d z 2 or- 
bitals pointed out in XAS and EELS experiments which 
shows 2-10% occupancy of d z 2 orbitals^2^S and 15% dop- 
ing dependent occupancy of p z orbitals^ in all hole doped 
High-T c compounds). Henceforth the multiband p — d 
model will be used. 

Let us consider the Hamiltonian with the following 
general structure: 

H = E (ex - rfnfxa + E E r #' c k^'* 

f,X,cr f^gX,X',cr 

+ \ E E vh x ' c )x ai cfx* 3 ci xla2 c g xi ai ,m 

f,g,X,X' ci j2 ,3,4 
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TABLE II: Hopping parameters and single electron energies 
for NCI2CUO4 obtained in WF projection procedure for differ- 
ent sets of trial orbitals (all values in eV). Here x 2 , z 2 , p x , 
p y , p z denote Cu-d x 2_ y 2, Cu-d 3z 2_ r 2, Ol-p^, 01-p y , 02-p z 
orbitals correspondingly. The 3-d and 4-th columns corre- 
spond to bases of the 3-band and the multiband p — d models 
respectively. 



Hopping 


Connecting 


Cu-2; 2 , 


Cu-x 2 , z 2 , 




vector 


O P X , Py 


O-Px, Py, Pz 






E x 2 =-1.989 


E x 2 =-1.991 






E Pa =-3.409 


E.2 =-2.778 






E Py =-3.409 


E Px =-3.368 
E Pi =-2.30 


', — PTT 

t(x 2 ,x 2 ) 


(1. 0) 


0.01 


0.01 


t'(x 2 ,x 2 ) 


(1. 1) 


-0.00 


-0.00 


t{z 2 ,z 2 ) 


(1, 0) 




0.01 


t'(z 2 ,z 2 ) 


(1 1) 




0.00 


t(x 2 ,p x ) 


(0.5,0) 


1.18 


1.18 


t'(x 2 ,p x ) 


(0.5,1) 


-0.06 


-0.06 


t'(x 2 ,Px) 


(1.5,0) 


0.04 


0.04 


t"'(x 2 ,p x ) 


(1.5,1) 


0.00 


0.00 


t(z 2 ,Px) 


(0.5,0) 




-0.29 


t'(z 2 , Px ) 


(0.5, 1) 




0.01 


t(2 2 ,Pz) 


(0, 0.5, 0.771) 




0.10 


t'(z 2 ,Pz) 


(1, 0.5, 0.771) 




0.02 


t(j>x,Py) 


(0.5,0.5) 


0.69 


0.67 


t'(Px,Py) 


(1.5,0.5) 


0.00 


0.00 


t(px ,Pz) 


(0.5, 0.5, 0.771) 




0.02 


t'(Px,Pz) 


(0.5,-0.5, 0.771) 




0.02 



part of the total Hamiltonian is diagonalizcd, allowing to 
classify all possible effective quasiparticle excitations in 
Cu02-plane according to a symmetry. To describe this 
process the Hubbard X-operators^ X™ <-> X p ^ q = \p) (q\ 
are introduced. Index m <-> (p, q) enumerates quasipar- 
ticle with energy oj m = e p (N + 1) — e q (N), where e p is 
the p-th energy level of the iV-electron system. There is 
a correspondence between Hubbard operators and single- 
electron creation and annihilation operators: 

c/ A CT =5>^( m ) x r> (11) 

m 

where 7a ct (to) determines the partial weight of a quasi- 
particle m with spin a and orbital index A. Using this 
correspondence we rewrite the Hamiltonian (|1L)|I 

H = ]T & - Nfx) XY + J2 E tjfxfxf. (12) 

f,V f¥=9rn,m' 

This Hamiltonian, actually, have the form of the multi- 
band Hubbard model. 

Diagonalization of the Hamiltonian i|l(J|) mentioned 
above gives energies e p and the basis of Hubbard op- 
erators XJ l . Values of the hoppings, 

t W = J- Tft'rt, (m) 7v<7 (m') , (13) 

a,X,X' 



where Cf\ a is the annihilation operator in Wannicr rep- 
resentation of the hole at site / at orbital A with spin a, 

n /Aa = c\ x<T C fXa . 

In particular case of cuprates and corresponding multi- 
band p—d model, / runs through copper and oxygen sites, 
index A run through d x 2_ y 2 = d x 2 and d 3z 2_ r 2 = d z 2 
orbitals on copper, p x and p y atomic orbitals on the Op- 
oxygen sites and p z orbital on the apical O a -oxygen; e\ 
- single-electron energy of the atomic orbital A. 
includes matrix elements of hoppings between copper 
and oxygen (t pd for hopping d x 2 <-> p x ,p y ; t pd /V3 for 
d z 2 «-> p Xl p y ; t' pd for d z 2 <-> p z ) and between oxygen 
and oxygen (t pp for hopping p x «-> p y ; t' pp for hopping 

Px,Py *-* Pz)- The Coulomb matrix elements Vf g in- 
cludes intraatomic Hubbard repulsions of two holes with 
opposite spins on one copper and oxygen orbital (U d , U p ), 
between different orbitals of copper and oxygen (V^, V p ), 
Hund exchange on copper and oxygen (J d , J p ) and the 
nearest- neighbor copper-oxygen Coulomb repulsion V pd . 

GTB mcthodiiSiSi consist of exact diagonalization 
of intracell part of the multiband Hamiltonian 11U|) 
and perturbative account of the intercell part. For 
La2- 2; Sr 2 ;Cu04 and Nd2- 2; Ce :r Cu04 the unit cells are 
CuC>6 and CUO4 clusters, respectively, and a problem 
of nonorthogonality of the molecular orbitals of adjacent 
cells is solved by an explicit fashion using the diagonaliza- 
tion in k-space^i. In a new symmetric basis the intracell 



are calculated straightforwardly using the exact diago- 
nalization of the intracell part of the Hamiltonian (|10f> . 

Again, in particular case of multiband p — d model, 
the essential for cuprates multielectron configurations are 
d w p 6 (vacuum state |0) in a hole representation), single- 
hole configurations d 9 p 6 , d 10 p 5 , and two- hole configura- 
tions d 8 p 6 , d?p 5 , d 10 p 4 , d 10 p 5 p 5 . In the single-hole sec- 
tor of the Hilbert space the b\ g molecular orbital, that 
we will denote later as \a) = {|f) , ||)}, has the min- 
imal energy. In the two-hole sector the lowest energy 
states are singlet state \S) with 1 A\ g symmetry, that 
includes Zhang- Rice singlet among other local singlets, 
and triplet states \T) = \TM) (M = +1,0,-1) with 
3 Bi g symmetryiSiSiii All these states form the basis of 
the Hamiltonian i|12|) . and they are shown together with 
quasiparticle excitations between them in the Fig. [S] 

In this basis relations (|llf> between annihilation- 
creation operators Cf\ a and Hubbard X-operators XJ L 
are 

c fPb<? = vX °f' T + 2<77 b Xj' , 

I /nvS-.TO v a,T2<y> 

C fPa<? = la(vv2X f -X f ' ), 
c fd z2 <? = 74 (J V2A / ' -X f ' ), 

c fp*<? = 7 P (cr\/2X ,' -X/ ), 
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|0> ± 



as/ / 

I \a 2 
T) 

! * I 7- 



|T) 
I S) 




FIG. 5: Schematic picture of states and quasiparticle exci- 
tations between them in Hubbard-type model 11121 . Here nh 
stands for number of holes, at numerates Fermi-type quasi- 
particles, states |0), \cr), \S), \T) represents basis oftheHamil- 
tonian 112C Also, bases of effective models are shown. 



and the explicit form of the Hamiltonian (|12fl is given by 

H P d = J2 [ £ i E x ? a + e ^xf 3 + x . 



TM.TM 
I 



A I 
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t f9 X f A 9 



J.SS Y S,S Y a,S 

l fg A f A g 



+2at% (Xf 
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TABLE III: Hopping parameters and single-electron ener- 
gies of holes obtained by fitting GTB band structure to ex- 
perimental data and in the ab initio calculations for p- and 
n-type cuprates (all values in eV). 







p-type 




n-type 




fitted 


ab initio 


fitted 


ab initio 
















£p* 


1.5 


0.91 


1.4 


1.38 




0.2 


0.14 


0.5 


0.79 




0.45 


-0.26 


0.45 


0.31 




1 


1 


1 


1 


tpp 


0.46 


0.63 


0.56 


0.57 


f> 

L pd 


0.58 


0.57 





0.08 


t' 


0.42 


0.29 


0.1 


0.02 



J' 



Xj 2 ^)(aV2X*' T0 - X a g < T2 °) 



Here a = —a. The relation between effective hoppings 
(|13fl in this Hamiltonian and microscopic parameters of 
the multiband p — d model is as follows^&ii: 



+oo _ 

t fg - 

j-SS _ 

l fg ~ 

t fg ~ 

,TT _ 

l fg ~ 



,ST 
r f<J 



-2tp d [j,fg2uv - 2t pp v fg v , 
-2t p dHfg{ v lx + ujb) - 2t pp Vf g v~f b , 



(15) 



2£ 

—?=-\ fg 2 lalz + 2t pp v fgl 2 a - 2t' pp \ fg 2 lplai 

^£fg"fz + 2t pp Xfg7a ~ 2t' pp t fglp . 



of annihilation-creation operators in the Hubbard X- 
operators representation^. 

Calculations^^ of the quasiparticle dispersion and 
spectral intensities in the framework of the multiband 
p — d model by the GTB method are in very good 
agreement to the ARPES data on insulating compound 
Sr2Cu02Ci2i2^i Other significant results of this method 

i) pinning of Fermi level in LSCO at low concentra- 
tions was obtained in agreement with experiments 42^°. 
This pinning appears due to the in-gap state, spectral 
weight of this state is proportional to doping concentra- 
tion x and when Fermi level comes to this in-gap band 
then Fermi level "pins" there. The localized in-gap state 
exist in NCCO also for the same reason as in LSCO, but 
its energy is determined by the extrcmum of the band at 
(tt/2, 7r/2) point and it appears to be above the bottom 
of the conductivity band. Thus, the first doped electron 
goes into the band state at the (w, 0) and the chemical 
potential \i for the very small concentration merges into 
the band. At higher x it meets the in-gap state with a 

c \ pinning at 0.08 < x < 0.18 and then /x again moves into 
J the band. The dependence fj,(x) for NCCO is quite asym- 
(14)metrical to the LSCO and also agrees with experimental 
datajSfi 

ii) experimentally observed^ evolution of Fermi sur- 
face with doping from hole- type (centered at (jt,tt)) 
in the underdoped region to electron-type (centered at 
(0,0)) in the overdoped region is qualitatively repro- 
duced; 

iii) pseudogap feature for LSCO is obtained as a lower- 
ing of density of states between the in-gap state and the 
states at the top of the valence band. 

In all these calculations the set of the microscopic 
model parameters, obtained by fitting to experimen- 
tal ARPES data ; 45 i 52 was used. Hoppings and single- 
electron energies are listed in Table IIIII values of 
Coulomb parameters are as follows: 



The factors fj,, v, A, £, x are the coefficients of Wan- 
nier transformation made in the GTB method and 
u, v, 7 X , 76, 7 a , 7 P , 7 2 arc the matrix elements 



U d = V d = 9, U p = V p = A 
J d = 1, J v = 0, V pd = 1.5. 



(16) 
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All results above were obtained treating the inter- 
cell hopping in the Hubbard-I approximation^ But the 
GTB method is not restricted to such a crude approx- 
imation. The Fourier transform of the two-time re- 
tarded Green function in energy representation can be 
rewritten in terms of matrix Green function D™™ (E) = 



CkXo 



'k\c 



1Xa (m) 7 L KW™ (E) 



The diagram technique for Hubbard X-operators is 
develop edS^4 and the generalized Dyson equation^ 
reads: 



Dka (E) 



G[°l (E) 



^ka{E)\ P k „{E). (17) 



Here, (E) and Pfc CT (E) are the self-energy and the 
strength operators, respectively. The presence of the 
strength operator is due to the redistribution of the spec- 
tral weight, that is an intrinsic feature of SCES. First 
time it was introduced in the spin diagram technique 
and called "a strength operator"^ because the value of 
Pka (E) determines an oscillator strength of excitations. 
It is also should be stressed, that T,ka (E) in Eq. I|17fl is 
the self-energy in X-opcrators representation and there- 
fore it is different from the self-energy entering Dyson 

equation for the Green function 



CfcAo 



J k\o 



Green function G k °^ (E) is defined by the formula 



4°J (E) 



= Go l (E)-P ka (E)t ko 



The 



(18) 



where G 1 (E) is the free propagator and ika is the inter- 
action matrix element (for the Hubbard model, t™™ = 
7ct (m) 7* (m') t k , and G? m ' (E) = 5 mm , / (iE - El )). 

In the Hubbard-I approximation at U ^ W the self- 
energy Efc CT (E) is equal to zero and the strength opera- 
tor P™ n (E) is replaced by P™« (E) -> P™« = 5 mn F™, 

where F™ = ^Xj?' p ^) + (^X'j' q ^ is the occupation factor 

So, in this approximation the following equation is de- 
rived from Eq. l(T7|) : 



D ( i}={G^{E)-P (T i ka } 1 P a 



(19) 



Using diagram technique for the X-operators it is pos- 
sible to find solution in the GTB method beyond the 
Hubbard-I approximation. But such discussion is far 
from the scope of this paper's goals. 

It should be stressed that the GTB bands are not free 
electron bands of the conventional band structure, these 
are the quasiparticle bands with the number of states in 
each particular band depending on the occupation num- 
ber of the initial and final multielectron configurations, 
and thus on the electron occupation. Bands with zero 
spectral weight or spectral weight proportional to dop- 
ing value x appear in the GTB approach. 



IV. LDA+GTB METHOD: RESULTS AND 
DISCUSSION 



In this Section we will describe the LDA+GTB method 
itself and some results of this approach. 

In LDA+GTB scheme all parameters of the multiband 
model are calculated within the ab initio LDA (by Wan- 
nier function projection technique, see Sec. Ill A|) and 
constrained LDA method^ Analysis of the LDA band 
structure gives the minimal model that should be used 
to describe the physics of system under consideration. 
Although LDA calculation does not give correct descrip- 
tion of the SCES band structure, it gives ab initio pa- 
rameters and reduced number of essential orbitals or the 
"minimal reliable model" . Then, the effects of strong 
electron correlations in the framework of this model with 
ab initio calculated parameters are explicitly taken into 
account within the GTB method and the quasiparticle 
band structure is derived. 

In Section [H] the ab initio calculations for undoped 
La2Cu04 and Nd2Cu04 are presented. One can see that 
in the 3-band model (Figs. 2] and |3J) it is possible to de- 
scribe the top of the valence band but not the lower ly- 
ing excitations withing 4 eV. The main effect of taking 
into account Cu-d 3z 2_ r .2 and O a -p z states for La2Cu04 
system is the proper description of the band structure 
(in comparison with LDA calculation) at the energies up 
to 4 eV below Fermi level (see Fig. 0). Of course, the 
ab initio LDA band structure is not correct in undoped 
cuprates, but it gives an indication what orbitals should 
be included in more appropriate calculations. Therefore 
if one needs to describe quantitatively the low-energy ex- 
citations of La2- z: Sr z: Cu04, the Cu-d 3z 2_ r 2 and O a -p z 
orbitals should be taken into account and the reliable 
minimal model is the multiband p— d model. In Nd2Cu04 
the Cu-ci 3z 2 _ r 2 and 02-p z states does not contribute sig- 
nificantly to the band structure (compare Figs. [3] and 0J) 
and the minimal model is the 3-band p — d model. Nev- 
ertheless to treat p- and n-type cuprates on equal footing 
later we will use the same multiband p — d model for both 
LSCO and NCCO with different material dependent pa- 
rameters. Hopping parameters decay rapidly with dis- 
tance (see Tables H] and |Hj so in GTB calculation we 
will use only nearest copper-oxygen and oxygen-oxygen 
hoppings which are listed in Table ITTT1 

In Refs. 15711581 ab initio calculations were done for 
YBa 2 Cu 3 07 and La2Cu04, and single-electron energy 
e Px = 0.9 eV was obtained. This value is very close to 
the one presented in Table IIIII But in Refsi 57 i 58 the Cu- 
s states were taken into account with energy e s = —6.5 
eV. Our LDA calculations shows that Cu-s bands con- 
tributes to the band structure shown in Figs. [21 and 0] 
at approximately 7 eV below and at 2 eV above Fermi 
level. Therefore Cu-s states does not contribute signif- 
icantly to the low-energy physics. But these states can 
contribute to the effective intraplane hopping parameters 
t' and t" between the nearest and next-nearest neighbor- 
ing unit cells. In our LDA+GTB method Cu-s states are 
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neglected. It could be a reason why for La2Cu04 our 
t'/t = -0.137 (see Table ITV) is less then t'/t = -0.17 
obtained in Ref. where influence of Cu-s orbital on 
hoppings was taken into account. 

There is a claim that pdn-honda^ and non-bonding 
oxygen states^ are very important in low-energy physics 
of High-T c cuprates. To discuss this topic lets start 
with analysis of ab initio calculations. Present LDA cal- 
culations show that anti-bonding bands 7r* of 7r-bonds 
(Cu-d+O-pn , see Figs. an d EJ) situated slightly be- 
low anti-bonding a a * bands of Cu-d 3z 2_ r 2+O a -p z origin 
in La 2 CuC>4 and slightly above anti-bonding a a * bands 
of Cu-<i3 Z 2_ r 2 origin in Nd 2 Cu04 (see Figs. El and 0J. 
GTB calculations^ show that states corresponding to 
a a * band contributes to the a\ g molecular orbital in the 
single-hole sector of the Hilbert space. This a\ g molec- 
ular orbital situated above b\ g state \a) = {| j) , | J,)} by 
an energy about 1.2 eV. From the relative position of 
a a * and 7r* bands in LDA calculations we conclude that 
the energy of molecular orbital corresponding to the w* 
band will be situated around energy of a\ g state. There- 
fore, it will be above \a) state by about 1.0-=- 1.4 eV. Also, 
both states corresponding to n* and <r a * are empty in un- 
doped compound and spectral weight of quasiparticle ex- 
citations to or from these states will be zero. Summariz- 
ing, 7r-bonds, as well as a a * states, will contribute to the 
GTB dispersion only upon doping and only in the depth 
in the valence band below 1 eV from the top. Moreover, 
since energy difference between triplet \T) and singlet 
\S) states is about 0.5 eV^ the contribution from the 
singlet-triplet excitations will be much more important 
to the low-energy physics. Although, 7r-bonds could be 
important for explanation of some optical and electron- 
energy loss spectroscopy experiments, but in description 
of low-energy physics of interest they could be neglected. 
The non-bonding oxygen states contribute to the valence 
band with energy about 2 3 eV below the top. That is 
why we will not take 7r-bonds and non-bonding oxygen 
states in our further consideration. 

Now we have an idea what model should be used and ab 
initio microscopic parameters of this model. As described 
in Section ITTT1 the GTB method is appropriate method 
for description of SCES in Mott-Hubbard type insulators 
and it's results are in good agreement with experimental 
data. Then it is natural to use this method to work with 
the ab initio derived multiband p — d model. 

The parameters <|15[) of the Hamiltonian in the GTB 
method derived from ab initio one are presented in Ta- 
bles \N\ and [3 for p- and n-type cuprates, respectively. 
Single-electron energies (in eV) and matrix elements of 
annihilation-creation operators in the X-operators repre- 
sentation were calculated for both LSCO: 

ei = -1.919, e 2S = -2.010, e 2T = -1.300, 

u = -0.707, v = -0.708, lx = -0.619, (20) 

7 6 = -0.987,7 a = -0.032, 7 P = -0.962, j z = -0.237, 



TABLE IV: Parameters of the multiband Hubbard model 
1121 and exchange integral J for LSCO obtained in the frame- 
work of the LDA+GTB method (all values in eV). Hoppings 
giving main contribution to the top of the valence band are 
shown by bold type. 
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(0,1) 


0.453 


0.679 


0.560 


0.004 


-0.086 
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(1,1) 
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0.001 


(0,2) 
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0.087 


0.002 


-0.016 


0.004 
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0.003 


-0.005 








-0.002 






TABLE V: The same as in TablelTVI but for NCCO. Hoppings 
giving main contribution to the bottom of the conductivity 
band are shown by bold type. 
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,00 

z p 


,ss 
z p 


,os 

i P 


,TT 
l P 


,ST 

l p 


Jp 


(0,1) 


0.410 


0.645 


-0.523 





-0.0052 


0.137 


(1,1) 


-0.013 


-0.076 


0.035 








0.001 


(0,2) 


0.058 


0.104 


-0.078 





-0.0002 


0.003 


(2,1) 


0.005 


-0.002 


-0.003 





-0.0004 






and NCCO: 

ei = -1.660, e 2S = -1.225, e 2T = -0.264, 

u = -0.756, v = -0.655, lx = 0.626, (21) 

75 = 0.984,7 a = -0.008, 7 P = 0.997, 7^ = 0.037. 

It is known 6 ! that sign of the hoppings in the t — t' — 
t" — J model changes during electron-hole transformation 
of the operators. Therefore, t p will have different signs 
in p- and n-type cuprates. In present paper we don't do 
electron-hole transformation of the operators and both 
t-t'-t"- J* and singlet-triplet t-t'-t" - J* models 
are written using hole operators. Because of that there 
is no difference in signs of the hoppings t p for the hole 
and electron doped systems presented in Tables HVl and 

m 

As the next step we calculate the band structure of 
the undoped antifcrromagnctic (AFM) insulating cuprate 
within the GTB method. Results for both GTB method 
with fitting parameters and LDA+GTB method with ab 
initio parameters ( Table ITTT^I are presented in the Fig- El 
for La 2 Cu0 4 and in the Fig. for Nd 2 Cu0 4 . 

The GTB band structure obtained for both phe- 
nomenological and ab initio sets of parameters is almost 
identical: the valence band, located below eV in fig- 
ures, and the conductivity band, located above +1.5 eV, 
divided by the insulator gap of the charge transfer ori- 
gin E ct w 2 eV; the undoped La 2 Cu04 and Nd 2 Cu04 
are insulators in both antifcrromagnctic and paramag- 
netic states. In-gap states at the top of the valence 
band and about the bottom of the conductivity band 
are shown by dashed lines. Their spectral weights and 
dispersions are proportional to doping x and concentra- 
tion of magnonaS. Therefore, for undoped compounds, 
in the Hubbard-I approximation used in GTB method, 
these states are dispersionless with zero spectral weight. 
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FIG. 6: The AFM band structure of La2Cu04 obtained in the GTB method with phenomenological set of parameters (left) 
and in the LDA+GTB method (right). In the left figure bands are labelled; in other GTB band structure figures relative 
positions of bands are the same. 
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FIG. 7: The AFM band structure of Nd2Cu04 obtained in the GTB method with phenomenological set of parameters (left) 
and in the LDA+GTB method (right). 



The valence band have bandwidth about 6 eV and con- 
sists of a set of very narrow subbands with the highest 
one at the top of the valence band - the so-called "Zhang- 
Rice singlet" subband. The dominant spectral weight in 
the singlet band stems from the oxygen p-states, while 
for the bottom of the empty conductivity band it is from 
d 2; 2_ 1 ,2-states of copper. Both methods give small, less 
then 0.5 eV, splitting between the x Ai g Zhang-Rice-type 
singlet band and 3 Bi g narrow triplet band located below 
the singlet band (e.g. in the Fig.[7|for Nd2Cu04 it is lo- 
cated at -1.5 eV). The energy of Cu-d 3z 2_ r 2 orbital plays 
the dominant role in this splitting in the GTB method. 
For La 2 Cu04 energy £^ 2 is smaller then for Nd2CuC>4 
(see Table ITTT1 . This results in smaller width of the sin- 
glet band for the LSCO compared to the NCCO: about 
0.5 eV and 1 eV correspondingly. 

However, for La2Cu04 minor discrepancies occurs in 
the dispersion of the bottom of the conductivity band 



near (n, 0) point obtained by GTB with phenomenolog- 
ical set of parameters and by LDA+GTB. This leads 
to the different character of the optical absorption edge 
in two presented methods. The absorption edge for the 
LDA+GTB is formed by the indirect transitions in con- 
trast to the GTB method with phenomenological set of 
parameters, where the momentum of excited quasipar- 
ticle is conserved by optical transition at the absorp- 
tion edge. For Nd 2 Cu0 4 both GTB method with fit- 
ting parameters and LDA+GTB result in the conductiv- 
ity band minima at the (7r,0) point (see Fig. Q). Also, 
in the LDA+GTB method the triplet band dispersion 
and the singlet-triplet hybridization are much smaller 
then in the GTB method with fitting parameters. This 
happens mainly due to the smaller values of t' used 
in LDA+GTB method, because it is this microscopic 
parameter that gives main numerical contribution (see 
Eqs. (HHJl, (HIJ) and J2TJ) to the iJJ and t s f ^ - hop- 
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pings that determines the triplet band dispersion and 
the singlet-triplet hybridization respectively. So, despite 
some minor discrepancies, both GTB method with phc- 
nomenological parameters and LDA+GTB method with- 
out free parameters gives similar band dispersion. 

Next topic that we will discuss in connection to the 
LDA+GTB method is the value of magnetic moment 
on copper Mcu- From the neutron diffraction studies 
of LaaCuO^S, and YBa 2 Cu 3 0g££ it is known that Af Cu 
is equal to 0.5/i^ where hb is Bohr magneton. There are 
two reasons of why Mq u is different from the free atomic 
value I.IA^b in 5* = 1/2 Cu 2+ , namely zero temperature 
quantum spin fluctuations and the covalent effect. Since 
each oxygen have two neighboring coppers belonging to 
different magnetic sublattices the total moment on oxy- 
gen is equal to zero. But due to p— d hybridization the p- 
states of oxygen are partially filled so these orbitals could 
carry non-zero magnetic moment Mo, while total mo- 
ment on oxygen will be equal to zero. Such space distri- 
bution of magnetic moment leads to the difference^ be- 
tween experimentally observed antiferromagnetic form- 
factor for La2Cu04 and the Heisenberg form- factor of 
Cu 2+ . In order to take into account covalent effects and 
zero quantum fluctuations on equal footing we will write 
down the expression for Mc u : 



M Cu = 2.28/! B (S z 



(22) 



where zero quantum spin fluctuations are contained in 
(S z ) and covalent effects are described by the weight u 2 
of the d 9 p 6 configuration. The last quantity is calcu- 
lated in the framework of the LDA+GTB method and 
equal to u 2 = 0.5. In paper |6^ the value (S z ) = 0.3 
was obtained self-consistently in the effective quasi-two- 
dimensional Heisenberg antiferromagnetic model for typ- 
ical in Ld2CuO\ ratio 10~ 5 of the interplane and intra- 
plane exchange parameters. Close value of (S z ) = 0.319 
was obtained in Ref. |63 where also the plaqucttc ring- 
exchange was considered in Heisenberg Hamiltonian. Us- 
ing Eq. H22|) and above values of u 2 and (S z ) we have cal- 
culated magnetic moment on copper Mq u = 0.4/Lts, that 
is close to the experimentally observed Mcu = 0.5/is. 

Summarizing this section, we can conclude that the 
proposed LDA+GTB scheme works quite well and could 
be used for quantitative description of the High-T c 
cupratcs. The LDA+GTB scheme also can be used for 
wide class of SCES - cupratcs, manganites, and other. 



EFFECTIVE LOW-ENERGY MODEL 



model results in the effective Hubbard and the t 
modeli 20 . 22 i 71 i 72 i 73 



J 



As the next step we will formulate the effective model 
for the multiband p — d model. Simplest way to do it is 
to neglect completely contribution of two-particle triplet 
state 3 B\g. Then there will be only one low-energy two- 
particle state - Zhang- Rice- type singlet 1 A\ g - and the 
effective model will be the usual t — J* model. But in 
the multiband p — d model the difference ex — es between 
energy of two-particle singlet and two-particle triplet de- 
pends strongly on various model parameters, particularly 
on distance of the apical oxygen from the planar oxygen, 
energy of the apical oxygen, difference between energy of 
d z 2 -orbitals and d x i -orbitals. For the realistic values of 
model parameters Et — £s is close to 0.5 eV^di contrary 
to the 3-band model with this value being about 2 cV. To 
take into account triplet states we will derive the effective 
Hamiltonian for multiband p—d model by exclusion of 
the intersubband hopping between low (LHB) and upper 
(UHB) Hubbard subbands. These subbands divided by 
the energy of charge-transfer gap E ct ss 2 eV (similar to 
U in the Hubbard model) and using perturbation theory, 
similar to Ref. 16^. with small parameter W/U wc can 
derive separate effective models for UHB and LHB. This 
procedure is schematically shown in Fig. [SJ And, as one 
can see, since the UHB and LHB in initial model l|14|l are 
formed by different quasiparticles (namely, cxq for LHB 
and ai, a-x, a% for UHB in Fig. [SJ, the effective models 
will be different for upper (valence band, hole doped) and 
lower (conductivity band, electron doped) subbands. 

We write the Hamiltonian in the form H = Hq + 
Hi, where the excitations via the charge transfer gap 
E ct are included in Hi. Then we define an operator 
H(e) = Hq + eHi and make the unitary transformation 
H(e) = exp (— ieS)H (e) exp (ieS). Vanishing linear in 
e component of H(e) gives the equation for matrix S: 



Hi 



Hq, S 



0. The effective Hamiltonian is ob- 



tained in second order in e and at e = 1 is given by 



H = Hn 



Hi, S 



For the multiband p — d model 

(|14fl in case of electron doping we obtain the usual t— J* 
model describing conductivity band: 



/+3 



l f n 9 > 



(23) 



When wc are interested in the low-energy physics 
(like e.g. superconductivity) it is useful to reduce the 
microscopic model to more simpler effective Hamilto- 
nian. For example, for the Hubbard model in the 
regime of strong correlations the effective model is the 
t — J* model (t — J model plus 3-centers correlated hop- 
pings H 3 ) obtained by exclusion of the intersubband 
hoppings perturbativclypi2i§2i2i Analysis of the 3-band 



here H$ contains three-centers interaction terms given 
by Eq. (| A. 1|> . Sf are spin operators and nj arc number 

of particles operators. The Jf g — 2 (t/f) /E ct is the 
exchange parameter. 

For p-type systems the effective Hamiltonian has the 
form of the singlet-triplet t — J* model describing valence 
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FIG. 8: (Color online) Quasiparticle dispersion and corre- 
sponding densities of states (DOS) in the t — J* model cal- 
culated for different number of taken into account coordina- 
tion spheres (c.s.). Chemical potential shown by the straight 
horizontal line was calculated self-consistently assuming 15% 
holes doping. 



band: 



H eff = H +H t + ]T H eff3 

+ J f3 [SfSs - \ n f n 9 



(24) 



Three-centers interaction terms -ff e //3 are given by 
Eq. (|A.2|) . Expressions for Hq and H t are as the follows: 



H 



= E 



ei *7 + £2sX s f ' s + s 2 tY X 



TM,TM 
f 



M 



t fg A f A s 
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X 



S,5 
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TO v a,T2a 
' X 9 



) + /i.e.] } 



The resulting Hamiltonian 1)24(1 is the two-band gener- 
alization of the t — J* model. Significant feature of ef- 
fective singlet-triplet model is the asymmetry for n- and 
p-type systems which is known experimentally. So, we 
can conclude that for n-type systems the usual t — J* 
model takes place while for p-type superconductors with 
complicated structure on the top of the valence band the 
singlet-triplet transitions are important. 

Contrary to the multiband p — d model's parameters 
that fall with distance rapidly, effective model parame- 
ters do not decrease so fast. This happens due to weak 
distance dependence of Wannier functions that determine 
coefficients Hf g , Vf g , A/ 9 , £/ g , Xfg^ which, in turn, deter- 
mine distance dependence of effective model parameters 

To demonstrate the importance of hoppings to far co- 
ordination spheres (c.s.) in the Fig. [SJ we present the 
dispersion and DOS in the t — J* model with parame- 
ters from Tabic IIVI The electron Green function 117|) 
has been calculated beyond the Hubbard-I approxima- 
tion by a decoupling of static correlation functions that 



r a a 
'S 



includes short-range magnetic order: (^Xj a X g 

Cf g . Here n p is the occupa- 



— 2°/ff) \ A / A g 



S f S g 



tion factors of the single-particle state, Cf g = 2 

arc static spin correlation functions which were self- 
consistently calculated from the spin Green's functions in 
the 2D t — J models As one can easily see from Fig. [SJ 
the dispersion with hoppings only to nearest neighbors (1 
c.s.) and to next-nearest neighbors (2 c.s., the so called 
t — t' — J* model) is quantitatively different around (0, 0) 
point and qualitatively different around (tt, 0) point from 
the dispersion with 3 c.s. (t — t' — t" — J* model) and 
more coordination spheres taken into account. 

Recent ARPES experiments^ show that the Fermi ve- 
locity vf = Ep/kp is nearly constant for wide range of 
p-type materials and doping independent within an ex- 
perimental error of 20%. We have calculated this quan- 
tity in the t—J* model with parameters from Tabic llVl in 
the approximation described above. In the doping range 
from x = 0.03 to x = 0.15 our calculations give very 
weak doping dependence of the Fermi velocity. Assum- 
ing the lattice constant equal to 4A we have vp varying 
from 1.6 eVA -1 to 2.0 eVA -1 . Taking into account ex- 
perimental error of 20% our results is very close to the 
experimental one. 



VI. CONCLUSION 

The approach developed here assumes the multiband 
Hamiltonian for the real crystal structure and its map- 
ping onto low-energy model. Parameters of the effective 
model 1(15(1 are obtained directly from ab initio multiband 
model parameters. The sets of parameters for the effec- 
tive models ((23(1 and ((24(1 are presented in Tables llVl and 
[3 for p- and n-type cuprates, correspondingly. 

The effective low-energy model appears to be the 
t-f-t"- J* model (JZU for Nd 2 Cu0 4 and the singlet- 
triplet t-t' -t" - J* model (EU for La 2 Cu0 4 . There 
is almost no difference in the band dispersion with ad- 
dition of numerically small hoppings to 4-th, 5-th, etc. 
neighbors. 

Summarizing, we have shown that the hybrid 
LDA+GTB method incorporate the ab initio calculated 
parameters of the multiband p — d model and the ade- 
quate treatment of strong electron correlations. 
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APPENDIX: EXPRESSIONS FOR 3-CENTERS 
CORRELATED HOPPINGS IN EFFECTIVE 
MODELS 

In the t — J* model (|23|l the 3-centers correlated hop- 
pings are given by: 



H 3 SS 



A / A m A ff _A / A m A S J 



Z(J \ A f A m^g + A / A m A . 
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(A.l) 



The three-centers interaction terms H e f f% in the effec- 
tive Hamiltonian (|24|l are much more complicated then 
in the t — J* model due to additional triplet and singlet- 
triplet contributions: 



tt L fm L mg S S 
H eff3 = p ^3 
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